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Abstract: The time- varying density of D-branes and anti-D-branes in an expanding 
universe is calculated. The D-brane anti-brane annihilation rate is shown to be too 
small to compete with the expansion rate of a FRW type universe and the branes 
over-close the universe. This brane problem is analogous to the old monopole prob- 
lem. Interestingly however, it is shown that small dimension D-branes annihilate more 
slowly than high dimension branes. Hence, an initially brany universe may be filled 
with only low dimension branes at late times. When combined with an appropriate 
late inflationary theory this leads to an attractive dynamical way to create a realistic 
braneworld scenario. 
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1. Introduction 

In a previous article we showed how D-branes may have been produced in a tachyonic 
phase transition at an energy near the string scale 0. While providing a well-defined 
non-perturbative method of generating branes, this mechanism suffers from the typical 
problem common to all topological defect theories - too many defects are produced. 

The mass of a D-brane is proportional to its volume and is inversely proportional to 
the closed string coupling g s . Thus non-compact branes at weak coupling are extremely 
heavy. Compact branes which wrap compact cycles of a spacetime M are less heavy, 
but are still much heavier than fundamental string states like electrons and heavy gauge 
bosons, etc. Thus, a gas of branes, whether compact or non-compact, is likely to cause 
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the energy density of the universe to exceed the critical energy density and turn an 
expanding universe into a contracting one 0. 

However, if g s — > 0, branes do not pose many cosmological problems because 
they then decouple from gravity in bulk. This is albeit that the mass m^p ~ V p /g s 
diverges. Gravity couples to D-branes through closed string emission and absorption. 
The "vertices" (in field theory parlance) for such interactions are generated by closed 
string vertex operators, V(k), which are proportional to g s . Thus as g s — > 0, the vertex 
operators "tend to" zero, closed string emission/absorption is suppressed, and gravity 
decouples from the branes. This is similar to the fact that the vertex for electron- 
photon interactions is proportional to the fine-structure constant a; and if a — > then 
electrons and photons no longer interact. Thus interestingly, heavy branes do not affect 
the gravitational dynamics of the universe in this context. 

This can be seen in a different manner. The perturbation Sgoo of a massive object 
on the surrounding spacetime in linearized gravity is given by the Newtonian potential, 

^oo = -^g^O if 0.-0 (1.1) 

because the d+ 1-dimensional Newton's constant varies as G^+i ~ gl and mo p ~ l/g s - 
However, for reasonable values of g s (which is related to the Yang- Mills coupling on 
the branes) branes will distort any cosmological spacetime. It is therefore imperative to 
dilute their densities via brane anti-brane annihilation, inflation, or something similar. 
In this paper we try to determine how brane anti-brane annihilation affects brane 
densities. 

If extra compact directions exist then the density of branes and anti-branes at late 
times will significantly depend on the process of decompatification of the large dimen- 
sions or compactification of the small dimensions, the initial topology of spacetime, and 
more generally the evolution of the compact directions. Our results can be summarized 
as follows. 
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Suppose that all of the compact directions are stabilised before brane production 
and that four non-compact directions are expanding. Then because of strong brane 
anti-brane interactions, branes filling the non-compact directions and separated only 
in the compact directions will quickly annihilate. Such branes will not feel the cosmic 
expansion which typically causes interactions to freeze out below some temperature. 
Thus branes with p > 3 which completely fill the non-compact directions and possibly 
wrap some of the compact directions will not pose any cosmological problems. Interac- 
tions between branes with p < 3 will be frozen out, and these branes may pose severe 
cosmological problems. Branes with p > 3 which do not fill the non-compact directions 
will pose problems as well. 

If all of the dimensions, compact and non-compact, are expanding at the time of 
brane production then branes of all dimensionalities will soon dominate the energy 
density of the universe and cause it to collapse. This is because brane interactions will 
quickly be frozen out and they will be unable to annihilate. If none of the dimensions 
have been decompactified, such that all the directions are compact but expanding as 
in the Brandenberger-Vafa model, then a variety of possibilites are possible. If the 
wrapped branes can generate a positive pressure and confine the compact directions 
via some sort of dilaton dynamics, then all of the wrapped branes may annihilate. This 
will ameliorate any brane problem but also prevent wrapped branes from allowing four 
compact directions to grow large. If the compact directions are not stabilized, then 
interactions will be frozen out and branes will be littered throughout spacetime. 

One interesting, but expected result from our calculations is that high dimensional 
branes annihilate more efficiently and drop out of equilibrium later. This is essentially 
because they interact more vigorously with other branes and matter in the bulk. This 
means that a Type IIB universe initally filled with branes may migrate to a universe 
filled with only a few D3 branes and many Dl strings. This may be very attractive 
from a braneworld point of view. 

The plan of this paper is as follows. First we consider the idealized case of non- 
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interacting branes in a static universe; this is essentially a reworking of the heuristic 
arguments of from a somewhat more algebraic point of view. Then we consider the 
case of interacting branes in a static universe. Next, we move onto the more realistic 
case of interacting branes in an expanding universe. Finally, we end with a catalog 
of how our results apply to scenarios with differing numbers of compact directions, 
expanding directions, wrapped directions, etc.. 

1.1 Non-interacting branes in Flat Space 

If we ignore long-range interactions and the cosmic expansion rate H' 1 , then the an- 
nihilation rate of a system of branes in a static spacetime M, should depend only on 
the distinguishing features of the branes such as their dimensionality (p + 1), and the 
geometry, topology and dimensionality of the spacetime M. 

Consider a non-coincident Dp brane and Dp' anti-brane. Label them as brane A 
and brane B and their worldvolumes as A and B. A and B will annihilate if they are 
coincident, parallel, and if p — p'. Anihilation occurs because tachyons arise which 
eventually condense to the minima of their potential. The minima for coincident anti- 
parallel branes corresponds to the disappearance / confinement of all open string degrees 
of freedom - i.e. annihilation |4]]. Tachyons also arise if two branes intersect at angles ||. 
However, the minima of the tachyon potential in that case correspond to changes in 
the geometry of the branes and not annihilation. 

A Dp brane and Dp' anti-brane will only meet each other if their spacetime paths 
intersect at some point in time - i.e. if their worldvolumes intersect. We have the 
following set theory identities 



A+B=A\B+B\A+AHB 




where 



A\B = A — An B 



4 



B\A = B - AHB. 



(1.3) 



This allows us to write 

dim(A + B) = dim(A - A n B) + dim(S - A n 5) + dim(A n B) 

= dim(A)+dim(S) -dim(AnS). (1.4) 

where dim(A D -B) is appropriately defined to ensure that the dim operator is linear. 1 
We now use the definitions 

dim(A + B) < D- dim(A) =p+ 1; dim(B) = j/ + 1 (1.5) 
where .D is the dimensionality of the spacetime. We then conclude that 

dim(AnB) > (p+l) + (p' + 1)-D (1.6) 

Thus, brane A and anti-brane S will surely meet at some point in time if (p + p') > 
D — 2. Equation |1.6| provides a lower bound on the dimensionality of the intersection, 
and the actual dimensionality of the intersection will usually be larger. For example, if 
(p + p') = D - 2, then dim(A nB)>0. The only way dim(A n B) = and A n B = 
can occur is if the worldlines of A and B are roughly parallel - an unlikely situation. 
Otherwise, they will intersect in a non-zero dimensional space. In the special case 
when p = p', and the Dp and Dp are parallel, any intersection of A and B will be p+ 1 
dimensional. This is crucial, because annihilation requires the brane and anti-brane to 
be coincident and that dim(v4 fl5) = p + 1. 

1 For example, if dim(A) = dim(_B) > dim(yl + B)/2 and the worldvolumes A and B are non- 
intersecting and parallel, then ( |1.4| ) would imply that dim(v4 n B) > which cannot be true. This is 
ameliorated by defining the intersection of the worldvolumes of two parallel objects to be intersecting 
in M n . This is reasonable because generic such worldvolumes will not be parallel and will intersect. 
Alternatively, one can work in projective space where for example, two parallel lines always intersect. 
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In string theory, D = 10. Thus, D3 branes will (almost surely) 2 not meet. But, D5 
and D7 branes will meet and annihilate. D3 branes and lower dimensional branes may 
meet. But the number of such configurations is a set of measure zero. For example, 
intuition tells us that two non-interacting low dimension branes can on rare occasions 
simply "run into each other." This will occur if brane A and brane B live in a smaller 
space W C M, which is a proper subset of M. Then effectively, D = dim(W). For 
example, if the worldlines of two DO branes lie in a 2D plane in a 10D dimensional 
space, then W = M 2 ; D = 2; and unless the worldlines are parallel, the branes will 
meet. The probability of two randomly moving objects lying in a proper subset of R 10 
is small. This would require that a subset of the transverse coordinates of A in M and 
B in M coincide exactly. The number of such coordinates is dim(M — W). Thus, the 
probability, P of lying in the subset is of order 



I N dim(M-W) 



which is vanishingly small. The actual probability will be larger, since A and B cannot 
be completely random surfaces. They must lie within their respective lightcones. 

However, if the transverse dimensions to W are compact with a compactification 
radius of L then 



I N dim(M-W) 



which may still be small, but does not vanish. Thus, compactification helps branes 



meet. But eq. |1.6| is still the principal criterion determining when non- interacting Dp 
and Dp' branes will meet. 

However, in reality branes are rarely non-interacting. Only parallel branes of the 
same dimension with the same orientation feel no attraction or repulsion to each other. 



2 In the measure theoretic sense. 
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An inter-brane potential exists for branes at angles, branes of different dimensions, and 
branes and anti-branes 0, ||, 0- 

1.2 Interacting Branes in Flat Space 

Recall that the charge of a brane relative to another brane depends on the relative 
orientation. The extreme case is when a brane is rotated by an angle ir relative to 
another brane. This corresponds to an anti-parallel or anti-brane. Just as particles 
and anti-particles feel a mutual attractive force, so do branes and anti-branes. Two 
branes with zero relative rotation angle are parallel branes. Thus, it is not surprising 
that the intermediate case of a rotation by < 6 < tt interpolates between no force 
and the large attractive force of brane anti-brane interaction; i.e. branes at angles feel 
forces. However, at a few special angles the force vanishes and the configuration is 
stable. 

Branes of different dimension also feel forces because open string excitations con- 
necting the branes are tachyonic. For example, a Dp is attracted to a D(p + 2). The 
force disappears once the Dp brane lies on the worldvolume of the D(p + 2) brane. A 
Dp is repulsed by a D(p + 6) brane, and unlike the Dp — D(p + 2) system does not flow 
toward a stable system 0. Dp — D(p + 4) brane configurations are supersymmetric 
and feel no forces || 

We will focus on the force between parallel branes and anti-branes of the same 
dimension. Only such systems can decrease the brane /anti-brane number densities in 
the universe. 

The potential $(r) between a Dp brane at a distance r away from a Dp is the 
analog of the gravitational potential, 0(r) in four dimensions 

0( r ) = _^ ^ $( r ) = cons t a nt (1.9) 

where is Newton's constant in D spacetime dimensions. In place of the mass, the 
potential contains the tension r p . The quantity Gr,r p is analogous to the effective grav- 
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ity of cosmic strings, Gfi, where pi is the string tension. It measures the gravitational 
strength of the branes. The potential energy V(r), of a brane and anti-brane is there- 
fore V(r) = rriDp^{r) in analogy with the purely gravitational case. Since k 2 = 8nGio, 
we can write V{r) as M 



V(r) = -f£g* = -£ (1.10) 
Here, (3 is a numerical factor characteristic of D-branes || H . 



= 7r -(9- P )/2 r( !_£) (L11) 

A brane and anti-brane in relative motion will have velocity and acceleration de- 



pendent potentials as well |T0|, For example, the largest velocity dependent piece 
is 

i) 2 h 2 ifi 

The first piece appears because of supersymmetry breaking and the second is 
present in susy and non-susy cases. Since both terms change the potential energy 
1.10|) by no more than a small numerical factor (of the same sign), we will ignore 



in 



them. We will also ignore acceleration terms. Note, we have left out the rest mass 
2m Dp in fllTTo!) . 



If some directions are compact then the brane inter-brane potential will depend on 
the average distance between the branes and anti-branes. If the distance is comparable 
to the size of the compact directions then the branes will notice the compactness of the 
extra-dimensional space. Hence, to calculate the effective brane anti-brane potential, 
the effects of an infinity of image charges will have to be added. This will soften the 
potential. However, the inter-brane separation is unlikely to be close to the compactifi- 
cation scale r c , as a gas of branes with a non-neglible density n, will have an inter-brane 
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seperation of n~^ d ~ p ^ which will usually be much smaller than r c - and this is the case 
that we will investigate. 

Because of long-range interactions the naive dimension-based arguments in the 
previous section do not hold for more general brane systems. Long rang forces lead 
to infinite scattering cross-sections. Thus branes and anti-branes will disappear due to 
annihilation if the universe remains static for a sufficiently long time and if annihila- 
tion doesn't require impossible topology change of any wrapped D-branes (impossible 
winding number transitions). 

1.3 Branes in an Expanding Universe 

We can calculate the brane and anti-brane abundances in an expanding universe as 



in 



T4 |13|j . First assume that the branes and anti-branes form dilute gases. This 
assumption allows us to ignore correlations between brane and anti-brane positions 
arising from the physics of a tachyonic phase transition which may have created them. 
Next, assume that the brane density equals the anti-brane density to ensure charge 
neutrality. The time evolution of the branes/anti-branes is then given by 

Off) 

— = -nH(d H -p H )-Tn 2 . (1.13) 

Here T characterizes the annihilation process; dn/dt is the time derivative of the density 
of branes n; and du is the number of expanding dimensions. The constant pn is the 
number of expanding directions which the branes fill. If p s is the number of static 
dimensions which the branes fill (non-compact or compact), then 



P = Ph +Ps 



;i.l4) 



The quantity (<i# — pn) represents the number of expanding dimensions in which the 
branes move. The first term in (|1.13|) represents the dilution effect of cosmic expansion. 
The second represents brane anti-brane annihilation. 
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The dn expanding dimensions may be compact or non-compact. Suppose that they 
are all non-compact and that the remaining d — du spatial dimensions are compact. 
Then branes living only in the compact directions will still be diluted because their 
inter-brane separations in the non-compact directions will grow. The only way branes 
in the compact directions will not feel the first term in ( |1.13| ) is if they are all coincident 
in the non-compact directions which is very unlikely, or they competely fill the extra- 
dimensions such that du = Ph- For example, if du = 3, then a gas of D3 branes all 
lying in the four non-compact directions and separated only in the compact directions 
will not suffer cosmic dilution. If the extra dimensions are themselves expanding or 
contracting, then branes will obviously feel additional dilution/contraction. 

Interestingly, according to ( |1.13| ) branes which fill some expanding directions such 
that ph 7^ 0, will feel less dilution than completely wrapped branes (pn = 0). Thus 
interactions between completely wrapped branes might be expected to freeze out more 
quickly than between incompletely wrapped branes (pn 7^ 0). 

The kind of branes that can fill the spacetime will be determined by the K-theory 
charges of the spacetime. However, since D-brane charges have yet to be classified for 
cosmological spaces, we will simply assume that they are the same as in flat space, i.e., 
for Type IIA/B theory all even/odd dimensional branes are possible, etc. The kind of 
cycles which branes can wrap and the number of such cycles are given by the Betti 
numbers of the compact extra dimensions. The Betti numbers largely dictate which 
brane configurations are possible. For example, a surface with 62 = &i = possesses 
no one-dimensional or two-dimensional cycles. Hence, D2/D1 branes will only be able 
to exist as contractable two-dimensional/one-dimensional loops on the surface. 

A more common compactification like K3 x T 2 possesses 2-cycles from the K3 and 
T 2 , two 1-cycles from the T 2 , 3-cycles from combining a 1-cycle of the T 2 and a 2-cycle 
of the K3, and 4-cycles by combining different 2-cycles of the K3 or by combining a 
2-cycle of the K3 with the 2-cycle of the torus. Finally, 5-cycles also arise by combining 
the entire 4-cycle of the K3 with a 1-cycle of the torus, and 6-cycles by combining the 
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entire K3 with the 2-cycle of the T 2 . Thus branes can wrap 1,2,3,4,5, or 6 dimensions. 

Since branes and anti-branes will annihilate only if they are parallel and share the 
same transverse space, we will assume that all the branes are parallel and possess the 
same transverse space V±. We can then take the dimension of n to be [n] ~ [L]~ d± . 

We parameterize the annihilation rate T as 



where m is some mass scale and A p is roughly constant. 

Although in this paper we will not worry about the functional dependence of A p , 
it will depend on the following. It will depend on the distribution of winding numbers 
of each cycle for each type of brane/anti-brane. This is because a brane and anti-brane 
can annihilate only if their winding numbers, (rii, n w ), around w compact cycles 
are the same. Collisions between branes can lead to inter-commutation of the branes 
and a change in their winding numbers. For example, a brane wrapping say cycle 7 
with winding number u> 7 = 1 and a brane wrapping say cycle f3 with winding number 
W/3 — 1 can collide and produce a single brane wrapping both cycles 7 and j3 with 
winding number {w 1 ,wp) = (1,1) if the intersection number of 7 and f3 is nonzero. 
Thus the intersection numbers tell us which winding number transitions are possible 
and A p will generally depend on them. 

For example for the K3xT 2 example, the intersection number *(A r • B s ) of a r- 
dimensional cycle A and a s-dimensional cycle B is one for two 1-cycles, 2b 2 (K3) for a 
2-cycle and 1-cycle. Many intersection numbers are non-zero. This means that branes 
wrapping the various cycles of the compactification can collide to produce branes with 
very different winding numbers. If some cycles are isolated, then only branes wrapping 
those cycles can interact with branes wrapping the same isolated cycles, reducing the 
volume accessible to such branes. This means that branes and anti-branes wrapping 
identical isolated cycles interact more easily. However, if the number of branes and anti- 




(1.15) 
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branes wrapping an isolated cycle are different, then due simply to topology, complete 
annihilation of the branes and anti-branes can not occur. 

In general, collisions between branes wrapping compact cycles and branes lying 
along non-compact directions will lead to branes wrapping compact cycles and lying 
along non-compact cycles M. Thus, spacetimes with both compact and non-compact 
directions are unlikely to have branes lying along only the compact or only along the 
non-compact directions. 

A p will also be affected by the presence of any points/locii at which the compact 
surfaces degenerate on the base manifold. At such points the volume of the cycles 
decreases. Such points may thus act as accumulation points for branes and anti-branes 
wishing to minimize their energy. 

For example, K3 possesses 24 singularities where a T 2 fiber degenerates over a base 
B of the K3. Thus, the two dimensions of the branes wrapping the T 2 become very 
small at the singularities. Presumably at such singular points it is much easier for a 
brane and anti-brane to annihilate, affecting A p . 

If some directions are fixed then the number of expanding dimensions satisfies 
dn < d. If some of the compact directions are expanding then the expansion rate 
will be affected by the wrapped branes. For example, the negative pressure of a gas 
of wrapped branes will increase the expansion rate. If the dilaton is allowed to freely 
evolve then the branes may actually inhibit and halt the expansion ||. However, this 
halting will only aid brane anti-brane annihilation if all of the directions are compact, or 
if the branes are coincident in the non-compact directions, or they fill the non-compact 
directions. In general, expanding non-compact directions will still dilute brane densities 
even if the compact directions stop expanding. 

Thus, we take our spacetime to be a FRW universe. We fix the dilaton. We ignore 
the expansionary effects of the brane density as it increases the expansion rate and 
freezes out brane interactions more quickly, making a FRW universe a conservative 
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choice. 3 Then the Hubble rate, H, of such a radiation dominated universe is: 



rp Ji(d H + l)/2 

H = -- = — — - (1.16) 

where C 2 is proportional to the effective number of spin degrees of freedom. Integrating 
(|Ll3l) , we find flTJ 



1 1 AC /m p \(^-i)/2 

7(T) ~ 7M + w+p#+(l-dff)/2 w / 

where /(T) = n/T dH ~ PH and is proportional to the comoving density. Here Tj is the 
initial temperature at which annihilation begins and /(Tj) is proportional to the initial 
comoving density. If w+pn + {^ — dff)/2 < then the annihilation is cut off for T <C 7*. 
If instead w + p# + (1 — du)j1 > as in ( |1.15| ) then for T Tj the number of branes 
/(T) becomes independent of the initial temperature T i and approaches the limit 



The annihilation rate will depend on brane anti-brane forces and brane interactions 
with massless bulk fields like the graviton, dilaton, etc. If the bulk is thermalized, 
numerous closed string fields will populate it. As the branes move toward each other, 
bulk particles will scatter off of them slowing them down. To scatter a brane by a 
large angle will require many bulk-particle scattering events, after which the brane's 
velocity will be randomized apart from a drift velocity in the direction of the attracting 
anti-brane. Let A be the distance the brane can travel before it is scattered by a large 
angle (mean free path). Also let r c be the distance between a stationary brane and 



Higher dimensional branes p = —(p/d)P r , 
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anti-brane for which capture is almost sure. Then, there are two important regimes: 
(1) the capture distance exceeds the mean free path, r c > A; and (2) r c < A. 

In the first case the branes will diffuse toward the anti-branes, and once the branes 
reach within r c of the anti-branes they will continue to be scattered such that all veloc- 
ities transverse to the direct path toward the anti-brane are randomized and effectively 
zeroed away (impact parameter = 0). The branes will thus collide directly with the 
anti-branes. If the drift velocity of a brane is sufficiently slow, the brane will be captured 
by an anti-brane in a bound state and eventually annihilation will occur. Otherwise 
the brane will simply pass through the anti-brane. 4 

Assuming collision occurs upon capture, the annihilation rate will be the rate of 
branes diffusing toward the anti-branes. The rate per unit density T, will be 



n 

where F is the number of branes approaching the anti-branes per unit time. As in the 
diffusion of charged particles through a plasma, this rate is 



where r is the mean free time and is the important quantity determining the rate of 
brane collisions. 

In the second case, once the brane is within r c of an anti-brane, its transverse 
velocities will not necessarily be randomized and the brane may streak past the anti- 
brane. The only way for capture to occur in this regime is for the brane to radiate 
away its kinetic energy. A moving brane drawn to an anti-brane will radiate via a 
Larmor-type formula. However, radiative capture is very weak and most branes and 
anti-branes will survive as we show later. 

4 We thank Ian Kogan and others for clarifications regarding solitons passing through one another. 






(1.20) 
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At high temperature the bulk particles will be very energetic and the mean capture 
length will exceed the mean free path (case 1). However, as the universe adiabatically 
expands, the temperature will drop until such a time that brane-bulk interactions are 
frozen out at a temperature Tf. After T — Tf, the mean free path will exceed the 
capture length and capture will occur only by radiative capture. Hence, below Tf, 
brane anti-brane annihilation will become very infrequent. 

2. Mean Free Path and Freezeout 

The capture distance r c can be obtained by the virial theorem, KE~-PE. Since KE ~ 
mv D , and PE~ —h 2 /r 7 ~ p , we find that 



However, D7 branes are not asymptotically flat and possess a logarithmically di- 
verging inter-brane potential: PE~ — lnr. Similarly, D9 branes have a quadratically 
diverging potential: PE~ — r 2 . Thus, the capture lengths for D7 branes and D9 branes 
are infinite. In any case, the notion of capture length for D9 branes is not a meaningful 
since the branes fill the whole space. 

We can obtain the mean free path, A, for a brane to be scattered by a large angle 
by bulk fields like gravitons, dilaton, etc., by writing A = Vd p t, where r is the mean 
time between collisions in which the brane is scattered by a large angle and vd p is the 
brane velocity. If the branes move with thermal velocities then vo p = ^/T/rri£) P and 
A = ^jT/m Dp T. 

Now, in a time r, the number of particles which are scattered is 



where a is the cross section for a particle (graviton, dilaton, etc,) to scatter off of a 
brane and v a is its average velocity. A particle can scatter two ways. Suppose the brane 





N s = v a ran, 



a 



(2.2) 
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fills the X 1 , ...,X P directions. Because the brane is an extended object, a particle in the 
transverse directions may directly strike the Dp brane. Since the brane has a width of 
order l/m s , unless the incident particle wavelength is smaller than l/m s it will reflect 
off the brane (the transmission coefficient is zero) |H |. For a wavelength less than l/m s , 



the particle will pass through the brane. Thus, for temperatures below the string scale 
the cross section will be a product of the geometrical size in the X , ...,X P directions, 
which is essentially the volume of the brane V p , and the cross-section to scatter off of 
it by streaking by it, cr s . Thus, 

a = V p a s for p < D — 2 (2.3) 

Note, if p = D — 2 (e.g. a D8 brane), a = V p . 

Now, we want the time necessary for scattering by large angles. A particle can be 
scattered by a large angle when the amount of energy AE, imparted to it by thermal 
collisions is of order the mass of the target object, in this case - the brane mass, m£> p . A 
gas of thermal particles will impart an energy AE ~ T to the target particle. Thus the 
target particle will be scattered by a large angle after approximately rriop/T collisions 
with lighter (massless) particles. Thus setting N s = m Dp /T |L2| 



r = — ^ (2.4) 
v a an a T 

We have v a ~ 1 and n a = aT D_1 since D — 1 is the number of spatial dimensions 
in which the light particles move. Thus 

_ 1 rn Dp 

Note that A grows with m,D P since the heavier the D-brane is, the harder it is to deflect 
by bulk scattering. 

We can now write down an expression for the temperature at which the mean free 
path equals the capture length, Tf. 
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Setting r c in eq EJ] equal to A in eq. 2.5, we find 



(2.6) 



where 



* (D-l/2)(7-p)-l V {D-l/2){J- V )-l (2 - 7) 
We parameterize the scattering cross section, a s as 



B rm\i 

°> = -ZJ-T (2- 



m^" 1 VT. 

where m is some mass scale. As section |5] shows for closed string scattering off of 
branes, q = 2 and m = m s . That calculation was done in the zero recoil limit. For 
compact branes, recoil will be more important and q may deviate from its value q = 2. 



Using (|2.8|), we find 




!//(!-«!/) , s 



1-9!/ 



(2.9) 



The exponent in (|27 

i±i^ = L f — L- + K — !— ( — + -) (2.io) 

l-qy 9.5-^-9 \7-p 2^ 9.5-9^7"? 2/ 1 v 
The important feature of ( |2.10| ) is that the exponent is strictly positive, (p < 7). 
Using q = 2 we find 

T <~[vJ (2 ' n) 

To make the dependence of Tf on p more explicit, we set V p ~ l v and calculate dTf/dp. 
Here, £ is some length scale, and the binary relation ~ is used because the different 
directions may have different sizes. 



-17- 



dT f \n£ ( r, 1 i\ / i \ 713(7^+2) 

Thus dTf /dp < 0. This implies that the higher the dimension of a brane (the smaller 
the co-dimension), the lower its freezeout temperature is, since V p grows with p and 
the exponent in ( |2.11| ) also grows with p. This is particularly true if the cycles the 



branes wrap, or non-compact dimensions the branes fill, are much larger than the string 
scale. Low dimension branes drop out of thermal contact with the bulk fields first. 
Higher dimensional branes continue to interact with bulk fields at lower temperatures 
and are thus captured by the same dimensional anti-branes at lower temperatures. 
Intuitively, this is plausible because the brane anti-brane potential grows with p, and 
hence the capture distance grows with p. Also, higher dimensional branes can more 
easily interact with bulk fields because their extended area is exponentially larger. This 
is essentially a selection rule which states that a universe populated with branes will get 
rid of higher dimensional branes more easily. We expect this, since higher dimensional 
branes can meet much more easily in the non-interacting case and presumably also 
in the interacting case. This implies that co-dimension dominates over mass in the 
annihilation process. One might have expected larger mass objects to drop out of 
equilibrium first. However, because of the extended size of the objects, this appears 
not to happen. For branes wrapping small compact directions with radii ~ y/a' \ then 
V p ~ 1, weakening the selection rule. But, in most cases the radii will be at least tens 
of string lengths. This is because brane formation is a stringy process occurring around 
an energy ~ m s (energy scale of tachyonic phase transition). Because, of the energy in 
the momentum modes of the strings and branes wrapping the compact dimensions, the 
cycles will typically be expanding. Thus by the time the branes have formed and start 
to look for one another, the compact universe will have substantially grown and the 
radii will be 3> \/a'. Furthermore, if we assume that annihilation occurs in the loitering 



phase of dilaton dynamics as in [15|, then radii ^> v a'. This is because loitering occurs 
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after the compact universe has grown substantially larger than its early string length 
size. 

^From ( 2.11 ) we can also infer that for fixed p that Tf grows with the number of 



static dimensions, p s , which a brane fills. If we take all such static dimensions to be 
compact, then p s is the number of dimensions a Dp wraps. If we write V p ~ £ p R p ~ Ps 
where R(t) is a cosmological expansion factor, then 

T (*) =^ ?±1 > (2.13) 

i \R p t p ) d Ps lp v ; 

This is expected as ( |1.13|) shows that for fixed p that Dp branes are frozen out more 

quickly as pu decreases or equivalently as p s = p — pn grows. 

Using equations |1.19| , |1.20| , |2^ we find that 



A p (m s \D-2 



-1 \ rp 



(2.14) 

m s x w / 

Since du — Ph < D — 1 and du < D — 1, we find that the condition for the number of 
Dp branes to be independent of the initial number: w + Ph + (1 — dn)/^ > 0, is always 
satisfied since w = D — 2. Because, Tf is essentially the lowest temperature at which 
our annihilation mechanism works, we find that at late times the number of Dp branes 
is 



D -2+p H + {l-d H )/2 fm,\^*- ^ T/ \ fl - 2 +^+d-^)/2 



aw A. r c UJ (2 ' 15) 

Note that Tf/m s < 1 and the part of the exponent in ( 2. 151 ) which is (Tf/m s ) PH 
decreases as pu grows. 

This means that the more expanding dimensions the branes fill, the smaller their 
final quantities are. If no compact directions exist or the branes live only in the non- 
compact directions such that p = Ph, then at freezeout the number of high dimensional 
branes is lower than the number of low dimensional branes for equal initial number 
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of branes. If compact directions exist and the branes wrap p s directions, then since 
Ph = P — p s , fewer high dimensional branes still remain. Thus, in addition to the fact 
that Tf is lower for higher dimensional branes, this further confirms that the number 
of higher dimensional branes is exponentially suppressed compared to the number of 
lower dimensional branes. 

For a fixed p, since pn oc — p s as a Dp brane wraps more directions its chance 
of being annihilated decreases. This is because according to ( |1.13| ) the more compact 
directions a gas of branes fills (the fewer compact directions it wraps) the more diluted 
it becomes. In the limit that the branes fill all of the expanding directions such that 
Ph = djj, Q2.15Q can be shown to take the form of f p {Tf) ~ l/T(Tf)t(Tf). 

Because f p does not vanish remnants remain. This is troublesome because once T 
falls below Tf annihilation can occur only by radiative capture, an inefficient process. 
For example, a brane may be attracted to an anti-brane. However, to form a bound 
state, the kinetic energy of both branes need to be radiated away. In the case of 
monopoles, it is known that almost all orbits leading to capture are nearly parabolic 
(eccentricity = 1) at the point of closest approach. Most of the radiation is given off 
during closest approach and the energy loss can change an initially hyperbolic orbit with 
eccentricity slightly larger than unity to a bound state elliptic orbit with eccentricity 
slighter smaller than unity. Once captured in this bound state, the brane then spirals 
in towards the anti-brane. At stringy distances, tachyonic open strings appear leading 
to brane anti-brane annihilation. 

Branes and anti-branes possess Ramond-Ramond charges. Thus, it is not surprising 



that accelerating branes and anti-branes radiate and possess a Larmor type formula [|T6 
We write the total power radiated by an accelerating brane/anti-brane as 



P= ^ P ,Dtt(D)h 2 v 2 (2.16) 
7 Pjj d describes how the power radiated per unit solid angle varies with p, and Q(D) 
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is the total solid angle in D spacetime dimensions; v is the acceleration which can be 
found from the trajectory. 

We will calculate the radiative capture rate for D6 branes wrapped on six compact 
directions and moving in four non-compact directions. Thus, from a four dimensional 
point of view, the D6 branes behave like point particles. In this case, D is effectively 
equal to four and p = 6. The general case of Dp branes moving in D — 1 dimensions 
is not much more interesting. This is because lower dimensional branes have an inter- 
brane potential which falls off faster and can less easily form bound states. Also, higher 
dimensional branes like D7 branes are expected to disappear quickly since the inter- 
brane potential of D7 — D7 brane anti-brane pairs does not drop off with distance, but 
grows logarithmically, V ~ — In r. 

Because we are interested in only the gross properties of the motion, we will ap- 
proximate the trajectory for a given energy E and angular momentum I by the same 
trajectory as traced out by a particle with the same energy and angular momentum in 
four dimensions. We also approximate the orbit during closest approach by an ellipse 
with eccentricity near unity. 

The radiation of kinetic energy leading to capture has been calculated in the case 



of monopole anti-monopole annihilation before by [13], [L7[ [Tj|, and the following calcu- 
lations are virtually identical. The amount of energy radiated in a frequency interval 
(u,u + du) is dE r = uN(u)du. Here N(u)dou is the number of quanta radiated with 
frequency near u. It can be parameterized as 



N(u)dw = ^M(K 2 1/3 (x) + K 2 2/3 (x))xdx (2.17) 
where Ki/ 3 ,K 2 /3 are Bessel functions and I is the orbital angular momentum, and 

l 3 u v 2 l 3 u 
X = W^i = WE~ k ( } 

The expectation value of the total energy radiated is then 
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(E r ) = / N(u)du = 76 ' 4 2/5 E k (2.19) 
Jo v 1 

where is the kinetic energy. All the kinetic energy is radiated away for angular 
momentum values less than a critical value, l c 



U = ^ ( ^ ) (220) 



V 



2 



1/5 



Capture will occur for I < l c . This gives a classical cross section of 



I x2 



a = n — (2.21) 
The thermally averaged interaction rate cross-section is therefore (using v ~ \JT /rri£>§ 

) 

r = W«l>«^-(^)* (2-22) 



m D6 



where 



ir=(2 4 7r 7 7 6 ,3) 1/5 ^ 4 oc^ (2.23) 

Note that the mass term is not the string mass, m s , but rather the mass of the D6 
brane, which is not exceedingly large because the D6 was taken to be wrapped on a 
compact surface. 

Since in this case pn = 0, dn = 3, and w = 9/10, we have w + pn + (1 — dn)/2 < 0. 
We know from the discussion around ( |1.17| , |1.18| ) that the annihilation is therefore 
cutoff and that the brane and anti-brane densities do not appreciably change from 
their values at T = Tf. Thus, the number of branes is not driven to zero and brane 
anti-brane annihilation doesn't resolve the brane problem. For the opposite case where 
a D3 brane fills the expanding four- dimensional spacetime and moves through the static 
six dimensional compact space, the annihilation rate T will be given by an expression 
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similar to fl2.22j ) and the exponent w is expected to be positive. Then w + pn + (1 — 
dn)/2 = 14/10 + cLh/2 will be positive definite because pn = du- This means that 
because the D3 branes fill the expanding space their annihilation is not cutoff and 
they continue to annihilate until they completely disappear. This is despite the fact 
that brane-bulk interactions are frozen out because the temperature decreases as the 
universe expands. 

3. Catalogue of results 

Below we list a somewhat pedantic catalogue of results - how our results apply to differ- 
ent spacetimes. The spacetimes differ in the number of compact, static and expanding 
directions. 

1. All directions non-compact and expanding: Brane interactions are frozen out and 
branes are leftover. This corresponds to the case where dynamical compactifica- 
tion of the extra dimensions occurs after brane production. 

2. All directions non-compact and static: Interactions do not have to fight cosmic 
expansion. Thus all of the branes disappear if the directions are static for a 
sufficiently long period. This may correspond to an initial state where all the 
dimensions are initially non-compact, static and some directions are later dy- 
namically compactified. 

3. All directions non-compact, d 1 ^ are expanding, d n s c are static; branes fill only static 
directions: The distance between branes and anti-branes expands. Hence interac- 
tions are frozen out and branes do not annihilate fast enough. (This corresponds 
to the case of p < d™ c , or p < 6 if d 1 ^ = 4). 

4. All directions non-compact, d"]f are expanding, d n s c are static; branes fill only the 
expanding directions: This corresponds to p = d 7 ^ or p = 3 for d 1 ^ = 4. As 
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the branes are separated only in the non-expanding compact directions, they do 
not feel the expansion and interact unhindered by the expansion. Thus if non- 
expanding directions are static for a sufficiently long period, virtually all of the 
D(dfj — 1) branes disappear. 

5. All directions non-compact, dJ^f are expanding, d n s c are static; branes fill pu ex- 
panding directions and fill p s static directions: As long as the branes move in 
some expanding directions they do not annihilate efficiently. This only happens 
if p H <d%. 

6. All directions compact and static: There is no expansion for interactions to fight. 
All the branes disappear if the compact directions are static for a sufficiently long 
period. 

7. All directions compact and expanding: Expansion will eventually kill off interac- 



tions and brane anti-brane remnants remain. However, if as in @, 15] the dilaton 
is allowed to vary, then via some dilaton dynamics the branes may be able to halt 
the expansion of the compact directions for a period sufficiently long for them to 
annihilate. 

8. All directions compact, d c H are expanding, d c s are static; branes wrap only static 
directions: The branes effectively behave like point particles moving in a non- 
compact spacetime. They do not disappear. 

9. All directions compact, d c H are expanding, d c s are static; branes fill only expanding 
directions: This is again the case of expanding branes moving in a static space. 
The branes disappear. 

10. All directions compact, d c H are expanding, d c s are static; branes wrap pu expand- 
ing directions and wrap p s static directions: The compact directions which the 
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branes wrap will act to only dimensionally reduce the branes which move in the 
expanding directions. This then reduces to scenario 7. 

11. There are d c compact directions, d nc non-compact directions. The compact di- 
rections are static, and the non- compact directions are expanding. Branes with 
p < d nc feel the expansion and do not disappear. For p > d nc then: (a) if the 
branes completely fill the expanding directions (e.g. non-compact D3 branes in a 
4+6 split of spacetime) then they disappear; (2) if the branes only partially fill the 
expanding directions and wrap p w directions such p — p w < d nc then the branes 
feel the expansion and create a remnant problem; (3) in the unlikely case that all 
of the branes are coincident in the non-compact directions the brane separation 
is oblivious to the expansion and the branes disappear. 

12. This is the most general case. There are d c compact directions, d nc non-compact 
directions. d c H compact and d]£ non-compact directions are expanding. d c s compact 
and d™ c non-compact directions are static. Branes move in the expanding non- 
compact directions if p < dJjj, and thus do not disappear. If p > dljj then 
branes still do not dissappear if they do not completely fill the non-compact 
space, i,e., if p — p w < d 7 ^. If at the same time the branes wrap some of the 
expanding directions, then either the wrapped branes must halt the expansion of 
the expanding compact directions or the brane remnants will remain. If they do 
do halt the such expansion then the brane density tends to zero. 

13. The case of some compact contracting directions Contracting directions do not 
shrink distances in the expanding directions. Thus contraction is unlikely to aid 
brane annihilation. However, this deserves further investigation. 

4. Conclusions 

D-branes are attractive brane-world candidates. However, a number of fundamental 
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questions must be answered before they become persuasive alternatives to the standard 
big bang model. Several of these are: why does there appear to be only one braneworld? 
Why is that braneworld four-dimensional and not say a D5 brane? Why are D3 
branes favored over higher or lower dimensional branes? If braneworlds were littered 
throughout spacetime, they would typically over-close the universe. Hence, either other 
braneworlds are very far away, diluted by something like inflation, or forced to be very 
rare because of some dynamical anomoly/tadpole cancellation mechanism, or their 
interactions with matter in the bulk is highly suppressed, or there is only one brane- 
world - ours. 

We showed that it is very difficult to dynamically produce an universe which mi- 
grates from a universe filled with many branes to a universe filled with only a few 
branes. Brane anti-brane interactions are cutoff by the expansion rate and annihilation 
after that point becomes very inefficient. Dynamical brane anti-brane annihilation in 
a non-compact universe is in fact more unlikely because branes and anti-branes will 
typically not be parallel. Hence, brane anti-brane collisions will typically lead not to 
annihilation but a complicated configuration of intersecting branes. 5 

However, in a universe with d c extra dimensions which have already been stabilised, 
we showed that branes may disappear if they fill all of expanding directions, p > 
<i^ c . For four non-compact directions, this means that Dp branes with p > 3 may 
disappear. However, if the branes do not wrap all of the expanding directions, they will 
remain. For example, even D7 branes may remain if they are separated in only one non- 
compact expanding direction. Universes in which the expanding directions are compact 
may become brane-less only if the wrapped branes conspire to halt the expansion. 
Otherwise, they suffer from the same problems that universes with expanding non- 
compact directions suffer from. 

We found one interesting result: higher dimensional branes tend to annihilate 

5 On a compact surface, branes tend to align themselves so as to minimize their energies so that 
non-intersecting branes wrapping the same cycles become parallel. Thus this may not be true. 
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fastest, leaving lower dimensional branes behind. This is intriguing because it implies 
that in a scheme where branes are dynamically generated as in a tachyonic phase transi- 
tion that the branes which survive and cause cosmological problems are low dimensional 
branes like D3 and Dl branes (in type IIB theory). Higher dimensional branes dis- 
appear. This at least begins to answer the question of why the universe may be filled 
with D3 branes as opposed to higher dimensional branes. This means that only one 
additional mechanism - inflation - may be needed to make the brane-world idea much 
more persuasive. Such inflation would dilute the D3 branes, making an observer on 
any D3 think that there is only one braneworld. Also, it would dilute any line (Dl 
branes) or point (DO branes) defects on the surfaces of the D3 branes. Such inflation 
would also dilute any small dimensional branes in the bulk as well. It would be inter- 
esting to see whether this same result holds once we also include interactions between 
Dp — Dp' branes, where p ^ p' '. In general terms this problem is very complex as there 
will not only be forces between branes of different dimensionalities, but branes oriented 
at angles, and intersecting configurations attracting other intersecting configurations. 
We did not include these interactions in the paper as they are irrelevent to the central 
question asked in this paper: do branes and anti-branes annihilate fast enough to make 
themselves small components of the energy density of the universe? 
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5. Appendix: Closed String Scattering off of D-branes 

Closed strings fields such as the dilaton, graviton, anti-symmetric tensor field, and 
massive modes can interact with a brane by exchanging closed strings. One can think 



-27- 



of this as spontaneously created open strings on the brane absorbing and then re- 
emitting closed strings and then subsequently annihilating themselves. This is the t 
channel point of view, see figure [I|. The dual s channel process is: a closed string strikes 
the brane, creates an open string excitation which then propagates along the brane and 



is eventually annihilated when it emits another closed string [19, 20 



/ 

/ 



We are concerned most with the t channel because, the t channel scattering ampli- 
tude blows up for small scattering angle. This divergence is also common to scattering 
by other solitons like monopoles. 

The form of the scattering amplitude can be 
guessed as follows. Because of the exchange of a 
massless particle, there will be a pole and A ~ 
1/t where t is the exchanged momentum squared. 
i>r The numerator of the amplitude should possess 

\^ some momentum dependent factors because the 

(() 

\ interactions between the bulk and the brane occur 

\ via quadratic derivative interactions. 

\ Because the amplitude is stringy, there should 

h 

also be a series of Regge poles, allowing for not 
Figure 1: t channel scattering by a only massless closed string exchange, but the ex- 
brane. change of massive closed strings as well. Thus we 

expect that 



A=^ + ^ + ^ + --- (5.1) 
t t+1 t+2 v; 

where we have taken the poles, in appropriate units, to be n G 1,2,3.... = Z. The 
factors are quadratic in the initial and final momenta, pi,P2- The amplitude for 
massless exchange ai/t, should be proportional the g s , since for g s — > the brane 
should decouple from the bulk and no scattering should occur. The simplest way 
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to achieve this is for a\ ~ k 2 t p ■ (function of pi,P2), since the amplitude should be 
proportional to Newton's constant k 2 ~ g 2 . 

Because momentum is conserved only parallel to the brane, (in the limit of small 
recoil), the only invariants are 



s = 2a'(p?) N = 2q/(^)|| 

t = -a'{pi + p 2 ) 2 = -2a px ■ p 2 (5.2) 

More formally, in the limit of no recoil the tree level stringy amplitude can be found 
by calculating the disk amplitude with insertions of two closed string vertex operators 
in the interior of the disk. 

A= I ^^(V 1 (z 1 ,z 1 )V 2 (z 2 ,z 2 )) (5.3) 

J VcKG 

V\ and V 2 can be decomposed into two right moving operators and two left moving 
operators. After fixing the positions of three of them one finds the familiar Veneziano 
scattering amplitude |19], |20 



A = ^nf^. l^i " t*2)5 p+1 (Pi + pa) (5.4) 

r(i + s + 1) 

where (sai — £02) is a kinematic factor and is equal to sTr(ei ■ 62) if the initial and 
final state polarisations, e\ and €2, are orthogonal to the brane. The delta function 
insures that momentum is conserved on the brane. When expanded in the t-channel 
(f<l,s>l), the amplitude ( |5.4|) takes the form of ( |5.1|) . The constant c\ is 

c i = ^7- (-}) ■ {o! 2 ) ~ g s (5.5) 

The factor on the left is the normalization of the disk amplitude. Note that because 
the open string coupling constant g Q is basically the Yang-Mills coupling constant, in 
dimensions other than four it is dimensionful, 1/g 2 ~ a' 2 r p . The second factor in ( |5.5|) 
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comes from the insertion of two closed string vertex operators. The third factor comes 
from evaluating two sets of propagators like (X(zi)X(z2)) = a' In \zi — 22 1. 6 

The amplitude can also be derived from an effective gravity picture [2~0|. The 
NS-NS sector of the low energy action of Type II theories is 



I NS = J d w x^ 
The source term for a D-brane is 



-l~R--(V<f>) 2 --H 2 e-^ 
2n 2 2 v V) 2 



(5.7) 



I s = J d w x^g[S^B, u + + Sf/v] 



(5. 



To leading order the source functions will be delta function sources at the positions 
of the branes. For example, the Fourier transform of the dilaton source is S$(k) = 
—T p (p — 3)/v8. If one expands the metric as g^ v = r/ M „ + 2nh^ vi we can read off the 
dilaton-dilaton-graviton vertex V/^ as follows 



2Kh^d»<Pd V <P -> 2K€ flvP t lp^(p 1 ) ( j ) (p 2 ) 



and find that the amplitude 



(5.9) 



A = iS^G^Vh^ (5.10) 

is the same as the massless part of ( |5.1| ) and ( |5.4| ) if we perform a field redefinition 
h^ v — > 2nh^ u and <p ~^ V2k,(J). 

The high energy limit, s, t > of the amplitude (|5.4j) is exponentially damped 
because of the Regge poles (.4. ~ sexp(— a'sm(<f>)/2) |^TJ. However, it peaks in the 
small momentum transfer limit i < 1, but large energy limit s ^> 1. In that case 



'Note the conventions, 



(W)V P 9s (2n) P V^ +1 ' 



K 2 = tMLfia/*. 



(5.6) 
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^=^ = 77^20 ( 5 - n ) 
t 2 sin | 

where we have assumed that momentum of the incoming particles to the brane such 
that s = a'E 2 . This is the massless term in our heuristic derivation (|5.1|) . The above 
amplitude is the zero recoil limit. In this case, p\ = P2, because as in usual Compton 
scattering for m Dp — > oo |2"2] 



m— >oo 

"^ i + ^a-cos^) ^" 1 (5 - 12) 

Thus t — pi ■ P2 — a'E 2 (l — cos 6). Because of this and momentum conservation, 
the final state |/) will be identical to the initial state The probability of scattering 
is 



ms\j)\ 2 = W 9_I 

2p x -2pi E 2 sin 2 6/2 



m - i) = tthPht^ = ~ ^ :ri „ (5-13) 



where we have explicitly included the normalization, = (pi\pi) = 2E pi . 

The cross section will be proportional to V{% — > i) once we average over the angle 
6. Thus 



°. ~ Ifc / -TWT, ~ £ (5-14) 



9lJ_ [*" d6 g 2 
£ 2 2vr J Bmin sin 2 6/2 ~ £ 2 

The dimensionful factors ignored in Q5.14|) are simply factors of a' and are not so 
important in determining the energy dependence of a s . The integral in ( |5.14| ) diverges 
for 9 m in = 0. Now a small deflection angle corresponds to a large impact parameter p3| . 



However, because the bulk fields scatter off other bulk fields as well as the brane, the 
maximum impact parameter of a particle which solely scatters off the brane is bounded. 
A bulk particle will feel a long range force from other bulk particles. The momentum 
change due to brane-bulk scattering will be (see figure H) 
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Ap = / Fdt 



dt r 00 dt 



r 8 - p Jo (b 2 + t 2 )( 8 -p)/ 2 
# 1 r d{t/b) 

~ W-V J ((t/b) 2 + 1)(8-P)/ 2 ~ W-P 1 ' 

Figure 2: Scattering of a bulk par- where 6 is the impact parameter and t is time. The 
tide by a brane. Note, t is time, and angle of deflection will be 9 = Ap/p. If the mean 
r is the distance from the brane to free path of the bulk fields is £, and p ~ T, then 
the bulk particle. Q min ~ 1/(T£ 7 ~ P ). Hence, the forward scattering 

limit of ( |5.14| ) is not singular. 
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